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ABSTRACT 

The rotation curves of spiral galaxies are one of the basic predictions of the cold dark 
matter paradigm, and their shape in the innermost regions has been hotly debated over 
the last decades. The present work shows that dark matter annihilation into electron- 
positron pairs may affect the observed rotation curve by a significant amount. We 
adopt a model-independent approach, where all the electrons and positrons are injected 
with the same initial energy i^o ~ TtdmC^ in the range from 1 McV to 1 TeV and the 
injection rate is constrained by INTEGRAL, Fermi, and HESS data. The pressure 
of the relativistic electron-positron gas is determined by solving the diffusion-loss 
equation, considering inverse Compton scattering, synchrotron radiation. Coulomb 
collisions, bremsstrahlung, and ionization. For values of the gas density and magnetic 
field that are representative of the Milky Way, it is estimated that pressure gradients 
are strong enough to balance gravity in the central parts if £^0 < 1 GeV. The exact 
value depends somewhat on the astrophysical parameters, and it changes dramatically 
with the slope of the dark matter density profile. For very steep slopes, as those 
expected from adiabatic contraction, the rotation curves of spiral galaxies would be 
affected on ^ kpc scales for most values of Eq. By comparing the predicted rotation 
curves with observations of dwarf and low surface brightness galaxies, we show that the 
pressure from dark matter annihilation may improve the agreement between theory 
and observations in some cases, but it also imposes severe constraints on the model 
parameters (most notably, the inner slope of halo density profile, as well as the mass 
and the annihilation cross-section of dark matter particles into electron-positron pairs) . 
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1 INTRODUCTION 

Although cosmological observations are providing increas- 
ingly convincing evidence that non-baryonic cold dark mat- 
ter (CDM) is the building block of structures in the Uni- 
verse, the exact nature of dark matter remains a mystery. 
A promising approach to the problem relies on the possi- 
bility that dark matter particles annihilate into observable 
products such as photons, neutrinos, protons, anti-protons, 
electrons, and positrons. Thus one can aim for indirect dark 
matter detection by l ooking for signatures of the annihila- 
tion products (see e.g. iBertone et al.ll2005l ). 

The most commonly studied signature is the 
emission of gamma rays from the Galactic cen- 
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tre (e.g. iBcrgstrom e t al.l Il998l: iBahz fc Edsio 'l999|; 
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Stoehr et ai] l2003l: [p7ada et all 12004 nCesarini et all 12004: 



Aharonian et al.ll2006l: ISpringel etlolbOQSl: ICirelli fc Paiidl 



20091 : iFornasa et all l2009l : iBernal fc Palomares-Rui j bOloT 
among many others) as well as photons at other fre- 
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As the annihilation products travel through the 
surrounding medium, they heat and ionize the gas, po- 
tentially lea ving an imprint on the cosmic microwave 
back gr ound |Chen fc Kamionkowskil 2004 : I Colafrancesco 



20041 : IPadmanab han fc Finkbeineij l200a iMapelh et al 
20061: [ziia ng ct ail I2OO6I. l2007l: JHooper et all l2007l: 



Cirelli et all l200ii: JGaUi et all l2009l : ICumberbatch et ahl 
2OO9I : iKanzaki et al.ll20ld ) and the HI 21 cm spectral line 
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l|Furlanetto et al.ll2006l : IValdes et al.ll2007l : IChuzhovl[20o3 ). 
Heating and ionization of the surrounding baryonic gas can 
also affect the formati on of the first s tars 1 focc o et al.ll2008l: 
Spolvar et al.l i2008: Nataraian et al.ll2009 : Ripam o nti et al 



201C ) and the formation and evolution of galaxies ([Ascasibar 



2007I - iRipamonti et al.ll2007l : iNataraian et al.ll200a i 



In this paper, we investigate the contribution of dark 
matter annihilation to the total gas pressure and consider 
the possibility that it has a significant effect on the rotation 
curve of spiral galaxies. While rotation curves provided one 
of the first and most important pieces of evidence fo r the 
existence of dark matter (see e.g. ISofue fc Rubinll200ll . and 
references therein), their shape in the inner regions of gas- 
rich dwarf and low surface brightness (LSB) spiral galaxies 
is o ne of the outs tanding issues in modern cosmology (see 
e.g. Ide Blokll2010l . for a recent discussion). 

Observationally, rotation curves are found to rise 
approximately linearly with radius, consistent with 
a constant de nsity core in the dark matter distri- 
bution (e.g. 



Flores fc Primack 
Kravtsov et al.l 
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Burkerd 1 19951: 

200ll; Ide Blok et al 

Marchosini ct al. 200 2.: iGentile et al. 

20041: ISpano et all I2OO8I: [Kuzio de Narav et all 120081 : 

de Blok et al.ll2008l : [Oh et al.ll2010l ) rather than the steep 



20041: iDonato et all 



pow er law predicted by cosmo l ogical N-body simu l ations 
(e.g. 'Dubinski & Carlbcre"l99l'; 'Navarro et al."l996', 'l997' 



2004 
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2010: Fukushige fc M akino 1997 
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Fukushige et al 
2OO5I : iDiemand et al 



2004; H avashi et al.l |2004|: iReed et al 
20051120081 : lAscasibar fc Gottlobei 



20081 : iGao et al. 2008 ) . Several modific ations to the stan- 
dard CDM scenario, suc h as warm JColin et al.l 2000l: 
Sommer-Larse n fc Dolgov* |200l|), repulsive JGoodmarJ 



2000l}. flu id (fPccblcs 2000), fu zzv dHu eraIll200Cll. de cay- 
ing (Con 2001), aimi hilating JK aplinghat et al. 2 000l'). or 
self-i n teracting (Spergel fc Steinhardt 2000 : Yoshida et al.l 
l200d : IPave et al.l I2001I ) dark matter, and even alternative 



theories of g ravitation (e.g 

ISanders fc Mc Gaugh 2002 : iGentile et al. 201 



McGa ugh fc d e Blokl 1 19981 : 
have been 



proposed in order to explain the discrepancy. 

Here we focus on the energy density associated to elec- 
trons and positrons arising from dark matter annihilations, 
neglecting other processes, such as dark matter decay, or 
other annihilation products, such as protons and antipro- 
tons (whose contribution i s severely constraine d by recent 
observational data; see e.g. lAdriani et al.ll2009br ). We adopt 
a model-independent approach, in which all particles are 
created with the same initial energy Eo ~ TTidmC^. Results 
for a particular dark matter candidate can be obtained by 
convolution with the appropriate source function. 

Since the characteristic energies involved are of the or- 
der of the mass of the dark matter particle, and this mass 
is usually much larger than the rest mass of the electron, 
electrons and positrons will be relativistic at the moment 
of their creation. However, they can efficiently lose their en- 
ergy through different processes, such as inverse Compton 
scattering (ICS), synchrotron radiation. Coulomb collisions, 
bremsstrahlung, and ionization. Throughout this paper, we 
will often use the Lorentz factor 7 to express the energy 



E = 77TicC^, where me denotes the rest mass of electron, 
and c is the speed of light. 

The pressure associated to these particles, hereafter re- 
ferred to as "dark matter pressure", is given by 
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where the electron-positron spectrum ^ (r, 7) is the number 
density of particles with Lorentz factor 7 at a radius r from 
the centre of the dark matter halo. The pressure gradient 
induces an acceleration 

dPdm(r) 



ddmir) 
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dr 
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Ps{r) 

where pg{r) is the gas density at radius r, that opposes the 
gravitational pull towards the centre, affecting observable 
quantities such as the circular velocity 
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It is our aim to show that, depending on the model pa- 
rameters, the contribution of dark matter pressure to the ro- 
tation curve may not be negligible. Section [2] describes the 
procedure followed to estimate the electron-positron spec- 
trum. The ensuing dark matter pressure is presented in Sec- 
tion [3l and the role of each astrophysical parameter (gas 
density and ionization fraction of the interstellar medium, 
intensity of the magnetic field, and dark matter density pro- 
file) is discussed in detail. The effect on the rotation curve 
is investigated in Section UJ and our main conclusions are 
briefly summarized in Section [5] 



2 THE ELECTRON-POSITRON SPECTRUM 

2.1 Propagation 

The propagation of electrons and positrons through the in- 
terstellar medium (ISM) is determined by the diffusion-loss 
equation 



d dn . . 
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0(x,7). 

We assume a diffusion coefficient of the form 
7^(7) = Ko-y' 
with Ko = 1.67 X lO^'^cm^s"^ and S 



(4) 



(5) 



0.7, independ ent of 

Galactic location (MED model in lDonato et al.ll2004l '). The 
energy loss rate 



&(x, 7) = - -^ (x, 1) = ^h (x, 7) 



(6) 



is a sum over the relevant physical processes, and the source 
term Q(x, 7) represents the instantaneous electron-positron 
injection rate. 

Given enough time, the electron-positron population 
will approach a steady-state distribution, ^ ^ (x, 7) — 0. 
Assuming that 6(x, 7) varies smoothly in space, the particle 
spectrum fulfills the relation 
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where 

2/(x,7) = 6(7) — (x, 7). 

Imposing ^ (x, 7) = at infinity, one obtains the Green's 
function 



(7) 



(8) 



G(x,7,Xs,7s) = 



exp 



/ |X-Xa|^ \ 
1^ 2AI5-J 



e(7 - 7e) 



(9) 



(27rAA2)^''^ 

and either the image charges method or an expansion over 
the eigenfunctions of the Unear differential operator may be 
used to derive t he Green's function for other boundary con - 
ditions (see e.g. iBaltz fc Edsio|[l999l : iDelahave et al.ll2009l ). 
The electron-positron spectrum is given by 
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where the quantity 

AA^ = a2(7) - a2(7s) 



d7s /d5c; 
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is related to the characteristic diffusion length of the elec- 
trons and positrons, and 7a denotes their initial energy. The 
variable A is defined as 



A' (7) = 



2J^(7) 
6(7) 



d7. 



(12) 



Considering the dark matter halo as a spherically-symmetric 
source, the spatial integral can be reduced to one dimension, 
and the electron-positron spectrum is finally given by the 
expression 

dn ^ 1 '^^p(-2i^) 

d7^'^'^^ 6(7) (27rr2AA2)i/2 



d7s / Ats Ts exp 



2AA2 



— exn 



''nAA^j"''p 



(-S)l«(-'^4^^3) 



Synchrotron radiation is another important loss mech- 
anism at high energies. The expression for the loss rate is 
similar to that of ICS, substituting the radiation energy den- 
sity in equation (|14|l by the energy density of the magnetic 
field, Ub ~ B^ /{Siv), where B is the intensity of the mag- 
netic field; 

4 (tt 



"syn(7) = |— ^7^%- 
6 rricC 



(16) 



For lower-energy electrons and positrons, Coulomb in- 
teractions with the thermal plasma m ust be taken in to ac- 
count. The loss rate is approximately (jRephaelill 19791 ) 



&Coui(7) ~ 1-2 X 10" 
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(17) 



where n^ is the number density of thermal electrons. 

Collisions with thermal ions and electrons also 
produce radiation through bremsstrahlung. The loss 
rate due to bremsstrahlu ng can be approximated as 
(|Blumenthal fc Gould|[l970h 



&brcm(7) ~ 1-51 X 10"'^ne7 [ln(7) + 0.36] 



(18) 



Additional energy losses come from the ionization of 
hydrogen atoms. The loss rate is given in lLongain (|l98ll ). 
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where nn is the number density of hydrogen atoms, q^ is 
the charge of electron, eo is the permittivity of free space 
and / is the ionization energy of the hydrogen atom. The 
number density of thermal electrons and neutral atoms can 
be expressed in terms of the total ISM gas density pg and 
the ionization fraction Xion as 



ne = ^Xio 



and 

WH = — (1-Xion) 



(20) 



(21) 



2.2 Loss rates 



respectively. 



Electrons and positrons can lose their energy by several 
physical processes as they move through the ISM. We con- 
sider ICS of cosmic microwave background (CMB) and 
starlight photons, synchrotron radiation. Coulomb colli- 
sions, bremsstrahlung, and ionization of neutral hydrogen 
atoms. 

The energy loss rates depend on the energy of the parti- 
cle. High-en ergy electrons and positrons mainly lose energy 
by ICS (e.g. ISarazinlll999l ). The relevant loss function is 



&ICS(7) = o 7^C^rad 

6 rricC 
where ctt is the Thomson cross section and 



C^rad = f/cMB + t/stars + Udu 



0.9 eV cm" 



(14) 



(15) 



is the combined radiation energy density of the 
CMB, starlight, and the rmal dust emission (see e.g. 
IPorter fc Moskalenkdl2005h . 



2.3 Source term 

Since the electrons and positrons in our model originate from 
the annihilation of dark matter particles, the production rate 
is dictated by the dark matter number density and the an- 
nihilation rate into electron-positron pairs, 

AN,± 



Q{r,^) = ndm(»-) ndm*('") (o-u). 



d7 



(22) 



where Wdm and nam* denote the number densities of dark 
matter particles and anti-particles, respectively, {av)^± is 
the thermal average of the annihilation cross-section times 
the dark matter relative velocity, and —^ — is the injection 
spectrum of electrons and positrons in the final state. 

Assuming all electrons and positrons are injected with 
the same energy 70 ~ mdm/mc, 



Q{r,l) = Qo{r) 5(7-70) 



(23) 
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2 X 10° 


1.022 X lO"-* 


2.6 X 10-30 


4.7 X 10-13 




2 X 10^ 


1.022 X 10^2 


2.6 X 10-28 


4.7 X 10-13 
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1.022 X 10-1 


2.6 X 10-26 


4.7 X 10-13 




2 X 10''* 


1.022 X 10" 


2.6 X 10-26 


4.7 X 10-21 




2 X 10'' 


1.022 X 10^ 


2.6 X 10-26 


4.7 X 10-23 
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2 X 105 


1.022 X 102 


1.0 X 10-25 


1.8 X 10-24 


0) 


2 X 10^ 


1.022 X lO''^ 


2.6 X 10-24 


4.7 X 10-25 


1 



Table 1. Initial Lorcntz factors, energies, cross-sections and pro- 
duction rates at the position of tiic Sun used in our calculations. 



where 
Qo{r) = 2 



Pdn-i{r) 



man 



(av) 



(24) 



is the local production rate per unit volume per unit time 
and ^(7 — 70) denotes a Dirac delta function. 

Although this is a rather coarse approximation, it has 
the advantage of being absolutely model-independent. More- 
over, the contribution of electrons and positrons to the gas 
pressure will be mostly determined by their total number 
and average initial energy, with the details of the injection 
spectrum playing only a minor role. The factor of 2 in equa- 
tion ((24} accounts for one electron and one positron pro- 
duced per annihilation event, and self-conjugate (Majorana) 
dark matter particles have been assumed. If dark matter 
particles and anti-particles were different, ridm = n-dm* ~ 
pdin{r) / {2mdin) and Qo would decrease by a factor of four. 
For the dark matter density pdmir), we consider a perfectly 
spherically-symmetric halo described by a density profile of 
the form 



Pdm(r-) 



i^n^+^j 



(25) 



where ra and ps denote a characteristic density and radius 
of the halo, respectively, and the a is the inner logarith- 
mic slope of the density profile. Local inhomogeneities that 
would boost the expected signal, such as small-scale dumpi- 
ness or the presence of subhaloes, are not taken into account. 
The shape of the dark matter density profile in the inner re- 
gions is far from being a settled question. As stated in the in- 
troduction, N-body simulations suggest that, at least in the 
absence of baryons, the profile should be quite steep near the 
centre (a ~ 1), in apparent contradiction with observations. 
Traditionally, it has been argued that the presence of gas and 
stars makes the pro file even steeper due to t he effects of adi- 
abatic contraction (jBlumenthal et al.lll986l ). although some 
rece nt claims have also been made in the opposite direction 
(e.g. lEl-Zant et al.|[20oil : iMashchenko et al'llSooel : lOh et al.l 
I2OI0I ). Given the current uncertainties, we have left the in- 
ner slope of the density profile as a free parameter of the 
model. 

For similar reasons, we also consider the injection en- 
ergy as a free parameter and investigate values of the initial 
Lorentz factor 70 between 2 and 2 x 10^, corresponding to 
energies Eq — ■yorricC^ from about 1 MeV to 1 TeV. The pro- 
duction rate Qo, on the other hand, is strongly constrained 
by different Galactic observations. At high energies, we con- 
sider observations of the electron-positron spectrum in the 
solar neighbourhood by HESS and the Large Area Telescope 
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Figure 1. Exclusion regions and production rates at the position 
of the Sun: black dots arc the adopted values of the instantaneous 
production rate Qg (see Tablc[TJ. Shaded regions above the black 
lines are excluded by INTEGRAL, Fermi and HESS data. 



(LAT) on board the Fermi satellite (JAharonian et al.ll2008l : 
lAbdo et al.ll2009l ). More specifically, the predicted amount 
of electrons and positrons cannot exceed the observed values 
for any Lorentz factor 7. Given the energy dependence of the 
observed spectrum, [^1 ~ E~^, and the energy losses, 
b{E) ^ E^, the most restrictive constraint comes from the 
spectrum near the injection energy, where propagation can 
be safely neglected and [^] « ^ oc E~'^. The max- 

imum production rate allowed by the data can then be ex- 
pressed as 



M^q) < fe(7o) 



dn 



dE 



(70) 

obs 



(26) 



Another, completely independent upper limit, valid at all en- 
ergies, can be obtained from the observed intensity of the 511 
keV line that measures the positron annihilation rate at the 
Galactic centre. In order to fully explain th e line with dark 
matter annihilations, it is necessary that (JAscasibar et al.l 
l200a) 



{c^v)s 



2.6 X 10-30 cm3 s-i 
so one just have 



1 MeVy 



Mre) < 2 



pdm{re) 



rridn 



Wv)s 



(27) 



(28) 



in order not to overproduce the observed signal. 

The corresponding exclusion regions are shown in Fig- 
ure [ij together with the production rates used in our cal- 
culation at the position of the Sun, r© = 8.5 kpc. These 
conditions constitute strict upper limits, since astrophysi- 
cal sources will also contribute to the relativistic particle 
budget, but the annihilation cross-sections they imply are 
comparable to or larger than the ones required to explain 
the cosmic dark matter density. 
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Hd 



(cru). 



(29) 



Therefore, we have set the injection rate according to 
equation \2ii) for m^^c? ^ 100 MeV, while a cross-section 
compatible with the relic density constraint, {av)^± — 2.6 x 
10~^^ cm'^ s~^, has been assumed for 100 Mev^ mdmC^ ^ 
10 GeV, and slightly larger values (based on the p ositron ex- 
cess observed by PAMELA; lAdriani et al.ll2009al ') have been 
used for rrtdmC^ ~ 100 GeV and 1 TeV. Numeric values are 
given in Table [T] 



3 DARK MATTER PRESSURE 

Apart from the initial energy and injection rate of the 
electron-positron pairs, related to the nature of the dark 
matter particle, there are many astrophysical parameters 
that determine the contribution of dark matter annihilation 
to the total gas pressure. We will first define a canonical 
model based on observations of the Milky Way and then 
investigate the effect of each individual component by vary- 
ing the values of the adopted parameters one by one. In 
all cases, we calculate the electron-positron spectrum as de- 
scribed in the previous section, and then estimate the dark 
matter pressure according to expression (ITJ. 

3.1 Canonical Milky Way model 

Our canonical mode l assumes a dark m atter density pro- 
file with a = 1 (JNavarro et all ligQ?! ). Vs = 17 kpc 
and ps = 6 X 10~^^ kg m"'', c onsistent with dynami- 
cal models of the Milky Way (e.g. iDehnen fc Binnevlll998l : 
iKlvpin et al.l 120021 ). The virial mass of the Galaxy is thus 
10^^ Mq, and the local dark matter density is pdm(?'Q) c^ = 
0.3 GeV cm~^. The ISM is mainly composed by neutral hy- 
drogen atoms (Xjon = 0) with number density pg/rrep ~ 
Icm"^ (|Dehnen fc Binnevl[l99i : iFerrier jlJOoH : iRobin et all 
120031 ). and it is permeated by a uniform, tangled magnetic 
field who se intensity i s B ~ 6 ^tG th roughout the Galaxy 
(|Ferriereli2001 ; Beck 200ll : lAscasibar fc Diaz 20ldi ]. 

The contribution of the individual loss terms described 
in Section 12.21 is plotted on the top panel in Figure [5] for a 
model similar to the canonical one, but with Xion = 0.5 
(in order to have a non-zero contribution from Coulomb 
collisions and bremsstrahlung). ICS and synchrotron emis- 
sion, being proportional to 7^, dominate at high energies, 
7 > lO''. Bremsstrahlung is important in the intermediate 
range 10^ < 7 < 10*, and Coulomb collisions and ioniza- 
tion, roughly independent on 7, dominate at low energies, 
7 < 10^ . The time taken by the electron-positron population 
to reach equilibrium is of the same order as the time 



treat (70) 



b(7) 



d7. 



(30) 



that the particles take to loose all their energy, shown 
on the middle panel in Figure [2] as a function of 70. Al- 
though this time may be larger than the orbital time at 
r ~ 500 pc (T ~ 10 Myr, assuming v ~ 220 km s^^) for 
JTidm > 100 MeV, a steady state will be reached as long 
as the conditions (dark matter and gas densities, magnetic 
field, etc.) evolve on timescales longer than ~ 100 Myr, 
which is relatively short in astrophysical terms. 
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Figure 2. Energy losses, equilibrium time scales, and charac- 
teristic diffusion lengths of electrons and positrons for pg/rrip = 



1 cm-'^', B = 6 fj.G, and Xi< 



: 0.5. On the top panel. Coulomb 



collisions, ionization, bremsstrahlung, synchrotron radiation, and 
inverse Compton scattering are represented by dotted, dot-long 
dashed, dot-short dashed, dashed, and solid lines, respectively. 



The steady-state electron-positron spectrum at the po- 
sition of the Sun, tq — 8.5 kpc, is shown on the left panel 
in Figure [3] for different values of the initial energy Eq. As 
stated above, all of the electrons and positrons are gener- 
ated with the same 70, according to Table [T] Propagation 
through the Galaxy and energy losses are accounted for by 
equation (|13p . The shape of the resulting spectrum is deter- 
mined by the value of 70, the production rate Qo{r), the loss 
rates 6(7) implied by the values of pg, B, and Xion, and the 
diffusion coefficient K{'y). Nevertheless, some insight may 
be gained by neglecting diffusion. For Kq — 0, 



dn 



6(7) 



(31) 



the electron-positron spectrum is almost flat when ioniza- 
tion dominates the energy losses, and there is a transition at 
7 ~ 10"^ (i? ~ 1 GeV) to the ICS-synchrotron regime, where 
^ ex 7"^. For low injection energies (1 MeV to 100 MeV), 
the normalizations of the spectra are identical because the 
value of Qo is only constrained by the INTEGRAL data, 
whereas other constraints impose lower values at higher en- 
ergies (see Figure [T}. In all cases, the spectra are cut at 
the injection energy 70 since no acceleration mechanism is 
included in our model. 

The electron-positron spectrum closer to the centre of 
the Galaxy (r = 10 pc) is shown on the right panel in Fig- 
ure[3l In general terms, the overall normalizations are higher 
than at the position of the Sun because of the higher dark 
matter density, and there is a sharp spectral feature near 
Eq. The characteristic diffusion scale AA plays an important 
role in both cases. As can be seen on the bottom panel of 
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Figure 3. Electron-positron spectrum at the position of the Sun (left) and at r = 10 pc (right) for our canonical model of the Milky 
Way and different values of the injection energy Bq- 
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Figure 4. Dark matter pressure as a function of radius in our canonical model compared with gas and magnetic pressure. 



Figure [2j A A depends on the Lorentz factor of the electrons 
and positrons. It is zero at the injection value, and it rapidly 
increases for lower energies until it saturates at a maximum 
value that depends on Eq. The spectrum at a given 7 probes 
the effective value of the production rate Q{r), averaged over 
the diffusion scale. This is not very relevant at the position 
of the Sun because the dark matter density does not vary 
much on kpc scales, but it becomes more important as one 
moves towards the central density cusp. For 7 <C 70, AA, 
and thus the effective production rate, is independent on 7. 
The larger Eo, the larger the smoothing scale, and there- 
fore the smaller the average density and the contrast with 
respect to the normalization at 8.5 kpc. As long as AA is 
constant, the shape of the spectrum remains the same, flat 



for low Lorentz factors and proportional to 7^^ in the in- 
verse Compton regime. Near the injection energy, AA be- 
comes very small, the effective production rate approaches 
the local source term Qo{r), much higher than the smoothed 
value, and the spectrum rises steeply just before the cutoff. 

Finally, the contribution of dark matter annihilations 
to the gas pressure at a given radius can be obtained by 
substituting the electron-positron spectrum in equation (|T]). 
Results for different values of 70 compared with gas and 
magnetic pressure in the galaxy are presented in Figure [l] 
compared to the thermal pressure of the gas Pgas = nkT 
(where n — pg/rup — 1 cm~ is the gas density, k is the 
Boltzmann constant, and we have assumed a temperature 
T = 100 K, appropriate for the neutral gas in the Galactic 
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Figure 5. Effect of varying the magnetic field intensity from 1 to 10 /iG (left), the ionization fraction X^^^ from to 1 (centre), and 
the ISM gas density from Pg/mp = 0.1 to 10 cm'^ (right). The pressure profiles for £0 = 1 McV, 100 McV, 10 GeV, and 1 TcV arc 
represented by dotted, dashed, dot-dashed, and solid lines, respectively. 



disk) and the magnetic pressure Pmag = B'^ /Sir, with B = 

For low injection energies, the pressure decreases 
sharply with distance from the Galactic centre. For Eo > 
1 GeV, diffusion keeps the electron-positron spectrum (and 
the ensuing pressure) roughly constant within a radius of a 
few kpc. The highest values of the dark matter pressure are 
found for an initial energy Eo — 100 MeV. Although the 
spectrum for Eo = 1 and 10 MeV is similar (approximately 
constant up to the cutoff at 70, because it dominated by ion- 
ization losses) and even has a higher normalization near the 
centre, the smaller upper limit of the integral in ((TJ yields 
significantly lower pressures. At high {Eq ^ 1 GeV) injec- 
tion energies, the dark matter pressure is also lower, due to 
the smaller number density of dark matter particles. Most 
of the dark matter pressure for an initial energy Eq, except 
Eo — 1 TeV, are higher than the pressure from the gas in 
the galaxy and for Eo — 10 MeV — 1 GeV, the dark mat- 
ter pressure is significantly higher than the pressure from 
magnetic fields. 



3.2 Astrophysical parameters 

We will now discuss the effect of the various astrophysi- 
cal parameters that enter our calculation of the dark mat- 
ter pressure, namely the intensity of the magnetic field, the 
density and ionization fraction of the ISM gas, and the in- 
ner slope of the dark matter density profile. As we did for 
the canonical model, we will compare the results obtained 
for different initial energies Eo from ~ 1 MeV to 1 TeV and 
vary each of the astrophysical parameters in turn in order 
to assess how much they infiuence the results. 

Magnetic fields affect the high-energy (7 > 10'^) tail of 
the electron-positron spectrum by setting the energy losses 
due to synchrotron radiation. As can be seen in Figure [2j 
in our canonical model with B — 6 fiG, the synchrotron 
term (|16|l is responsible for about 50 per cent of the energy 
loss at high energies, with ICS being responsible for most 
of the other 50 per cent. At low energies, energy losses are 
dominated by ionization of neutral hydrogen, and the con- 
tribution of synchrotron emission is negligible. The effect of 
varying B from 1 to 10 /iG is plotted on the left panel in 



PbC^ [GeV cm 



[kpc] 



1.0 


0.3483 


16.68 


1.2 


0.1975 


20.34 


1.4 


0.09946 


25.88 


1.5 


0.06658 


29.81 


1.7 


0.02469 


42.46 


1.9 


0.00615 


70.30 



Table 2. Characteristic density and radius of the dark matter 
density profile lj25|l as a function of its asymptotic logarithmic 
inner slope a. 



Figure [5l Not surprisingly, the results for an initial energy 
-Eo < 1 GeV are largely unaffected. At higher energies, the 
pressure at large radii decreases with the magnetic field in- 
tensity because of the more rapid energy losses. However, 
the diffusion length becomes shorter, increasing the effec- 
tive production rate and yielding a larger pressure near the 
centre of the Galaxy. 

As explained in Section 12.21 the density of the ISM gas 
and the ionization fraction A"ion regulate the energy losses 
by Coulomb interactions, bremsstrahlung, and ionization. 
For Xion = (our canonical model), the ISM gas is entirely 
composed of neutral hydrogen atoms, and the energy loss of 
the electrons and positrons with 7 < 10^ is dominated by 
the ionization process. At the other extreme, Xion = 1, the 
ISM is already fully ionized, and the relevant energy losses 
are Coulomb collisions and bremsstrahlung. Since the total 
loss by these processes is higher than the loss by ionization, 
the maximum pressure happens when Xion = (mid panel 
of Figure [SJ . The effect of changing the gas density from 
0.1 to 10 cm~^ is shown on the right panel. Higher densities 
yield lower dark matter pressures, simply because the energy 
losses are faster. 

Finally, we calculate the dependence of dark matter 
pressure on the inner logarithmic slope a of the dark matter 
density profile. When varying a, we also modify the charac- 
teristic density and radius in expression (|25p so that the dark 
matter density at the solar radius is equal to 0.3 GeV cm""^ 
and the virial mass of the Galaxy is 10^^ M©. The appro- 
priate values of pa and rs are quoted in Table [2] for several 
values of the inner slope. For a ^ 1.5, the production rate 
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Figure 7. Dark matter acceleration for different values of the injection energy Eq and the inner logarithmic slope a. 



Qo in equation (|24p diverges rapidly at r = 0, so we add a 
cutoff based on tfie local annihilation rate r(r), 



Qo{r) 



/5dm(j-)exp(-tor(r)) 



"Idn 



(av) 



(32) 



where to = f 3.7 Gyr is the age of the universe and r(r) — 

The dark matter pressure profiles obtained for a — 1, 
1.2, 1.4, 1.5, 1.7 and 1.9 are compared in Figure[6l Since the 
central dark matter density increases dramatically with the 
value of the inner logarithmic slope, this is, by far, the most 
relevant astrophysical parameter, only second in importance 
to the injection energy Eq related to the mass (and the pre- 
cise nature) of the dark matter particle. 



4 ROTATION CURVES 

The gradient of the dark matter pressure induces an accel- 
eration on the baryonic component that opposes the gravi- 



tational force. This acceleration, given by expression ©, is 
plotted in Figure [7] for all the injection energies considered 
in this work and compared to the gravitational accelera- 
tion g{r) = GM{r)/r'^ (represented by a thick solid black 
line). Each panel corresponds to a different value of the 
inner logarithmic slope of the dark matter density profile. 
For our canonical model with a = 1, g{r) ~ 27rGps'"s = 
1.3 X 10"^" m s~^ in the innermost regions, whereas in the 
general case described by equation ([25} gravity scales as 
g{r) oc r^~" for r <g r^. 

Depending on the model parameters, the acceleration 
caused by the electron-positron gas may be comparable to 
(or even higher than) the gravitational one in the central 
parts of the halo. For our canonical model (left panel of Fig- 
ure [7|, the pressure gradient is strong enough to overcome 
gravity for Eq < 1 GeV, and the radius at which both forces 
balance each other is of the order of 100 pc. The effect of 
dark matter annihilation is weaker, but perhaps still mea- 
surable, for Eq ^ 1 GeV. It would be extremely difficult to 
detect at 10 GeV, and completely negligible for larger parti- 
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8. Effect of dark matter annihilation on the rotation curve. From left to right, panels correspond to a = 1.0 (our canonical 
1.5, and 1.9, respectively. 



Name 


p^c'^ [GeV cm-3] 


r-s [kpc] 


XdBB 


v2 
Al00McV,HI 


2 
^IGcV.HI 


-^lOOMoV.HII 


v2 

-^IGoV.HII 


UGC 1230 


9.2799 


X 10-2 


15.152 


0.723 


0.713 


0.720 


49.825 


0.430 


UGC 5005 


8.5056 


X 10-5 


4206.9 


0.521 


0.520 


0.521 


4.726 


0.488 


LSBCF563-01 


4.5080 


X 10-^ 


188.61 


0.383 


0.381 


0.383 


27.368 


0.345 


UGC 4173 


8.5056 


X 10-5 


1765.8 


0.225 


0.225 


0.225 


1.223 


0.214 


UGC 3371 


8.5056 


X 10-5 


7302.7 


0.387 


0.383 


0.386 


17.775 


0.273 


NGC 1560 


8.5056 


X 10-5 


7131.5 


8.914 


9.183 


8.896 


659.459 


8.720 


DDO 189 


3.9944 


X 10-3 


135.71 


0.173 


0.171 


0.173 


40.124 


0.135 


NGC 4395 


3.8588 


X 10-1 


5.8478 


0.573 


1.627 


0.555 


743.952 


14.954 


NGC 3274 


4.4477 


X 10-1 


6.7189 


1.787 


7.139 


1.731 


2935.570 


33.516 


NGC 4455 


8.5056 


X 10-5 


8252.1 


1.428 


1.694 


1.419 


320.260 


1.643 


NGC 2366 


1.7431 


X 10-1 


2.3734 


1.246 


1.209 


1.245 


58.692 


1.204 


UGC 4325 


8.5056 


X 10-5 


23618 


1.326 


1.294 


1.300 


748.686 


4.459 


DDO 47 


8.5056 


X 10-5 


14107 


0.400 


0.303 


0.394 


260.841 


2.056 


DDO 185 


8.5056 


X 10-5 


5245.2 


2.036 


2.024 


2.026 


141.203 


1.896 



Table 3. Galax y nam e, characteristic dark matter density and radius, and reduced x^ values for the maximum disk models in 
Ide Blok &: Bosmal 1 120021 ). XdBB ''^^^1'^ to the original model without dark matter annihilation. The subscripts in the other columns 
denote the injection energy (100 MeV or 1 GeV) and the conditions in the baryonic medium (n = 1 cm— -^ and x^^-^ = for HI, 
n = 0.01 cm-'' and xion = 1 for HII). 



cle masses. These conclusions are very robust with respect to 
variations in the ionization fraction of the gas or the inten- 
sity of the magnetic field. The exact density of the interstel- 
lar medium has a somewhat larger influence on the results, 
partly because of its effect on the dark matter pressure (see 
Figure [5| and partly through the presence of the gas density 
in equation ([2]). For the extreme case Eq = 100 MeV and 
Pg/rrip = 0.1 cm-"', dark matter pressure is able to prevent 
gravitational collapse within the inner 2 kpc, compared to 
100 pc for a density of 10 cm-''. However, the qualitative 
picture is not changed. For Eq > 10 GeV, the gravitational 
acceleration dominates by several orders of magnitude at all 
radii, even for the most dilute gas. 

As shown in the previous section, the logarithmic slope 
of the density profile plays a critical role on the pressure 
profile. The accelerations for a = 1.5 and 1.9 are shown on 
the centre and right panels of Figure [T] respectively. When 
a — 1.9, the radius at which the relativistic pressure bal- 
ances gravity ranges from a few tens of parsec up to several 
kpc, and a sizable effect on the rotation curve of the Galaxy 
is expected for any value of the injection energy Eq < 1 TcV. 

Figure [8] shows the modified circular velocity profiles. 



according to expression ([3]), for different values of a and Eq. 
In our canonical model, the rotation curve of the Galaxy 
changes significantly for Eq = 100 MeV, it is slightly mod- 
ified for Eo = 10 MeV and almost imperceptibly for Eq = 
1 GeV. For higher values of the inner logarithmic slope, 
as those pred icted, for instance, in t he adiabatic contrac- 
tion scenario (|Blumenthal et al.lll98a ). it is more likely that 
the annihilation of dark matter particles leaves a clearly de- 
tectable imprint on the observed rotation curve. The scales 
on which such a signal would be visible, of the order of kpc 
in some cases, subtend several degrees on the sky for the 
Milky way, and may be observable as well in other nearby 
galaxies. 

In other to quantify the effect on the rotation curves 
of low surface brightness galaxies, we c ompute the model 
predic tions for the objects compiled by Ide Blok fc Bosmal 
(|2002h . Since it is not our aim to fit the data (which would 
require more careful modelling, beyond the scope of the 
present work), we simply take the observed rotation curves, 
as well as the quoted decomposition into stellar disk, gaseous 
disk, and dark matter halo. We consider their constant mass- 
to-light ratio and maximum disk models, adopting the cor- 
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Figure 9. Model predictions, compared to the observed rotation curves (data points with error bars). Solid lines are used for Eq = 
100 MeV and dashed lines for 1 GeV. Results for neutral and ionized gas are shown in black and blue colo ur, respectively. Orange lines 
depict the model without dark matter annihilation. The rotation curves by star and gas (data taken from Ide Blok fc Bosmall2002r are 
shown by red dot-dashed and black small dashed lines, respectively. 
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Name 


Psc2 [GeV cm-3] 


r-s [kpc] 


2 
>^dBB 


v2 
A100McV,HI 


v2 
AlGoV.HI 


-»^100MoV,HII 


2 
>^lGoV,HII 


UGC 1230 


3.7001 


X 10"^ 


9.0973 


1.091 


1.032 


1.072 


257.449 


5.447 


UGC 5005 


9.1650 


X 10'^ 


75.759 


0.175 


0.172 


0.174 


18.571 


0.102 


LSBCF563-01 


9.9553 


xlO-2 


20.850 


0.370 


0.355 


0.367 


155.857 


1.492 


UGC 4173 


8.5056 


X IQ-s 


2413.7 


0.124 


0.123 


0.124 


3.687 


0.110 


UGC 3371 


8.5056 


X IQ-s 


10073 


0.280 


0.273 


0.278 


36.263 


0.118 


NGC 1560 


5.6447 


X IQ-^ 


166.59 


2.330 


3.810 


2.315 


1905.660 


17.378 


DDO 189 


2.2759 


X IQ-^ 


7.4629 


0.109 


0.088 


0.105 


273.964 


5.120 


NGC 4395 


4.4477 


X IQ-i 


6.3057 


0.644 


2.476 


0.613 


1149.580 


30.584 


NGC 3274 


2.2174 


xlO-o 


2.9131 


0.941 


29.406 


0.808 


9205.720 


371.555 


NGC 4455 


8.5056 


X 10-5 


9702.7 


0.614 


1.078 


0.605 


466.972 


2.482 


NGC 2366 


3.8588 


X IQ-i 


4.1902 


1.935 


2.222 


1.900 


950.701 


16.488 


UGC 4325 


8.5056 


X 10-5 


36371 


1.096 


1.081 


1.047 


1737.360 


34.844 


DDO 47 


8.5056 


X 10-5 


16137 


0.272 


0.180 


0.265 


344.709 


3.186 


DDO 185 


8.5056 


X 10-5 


7876.7 


2.162 


2.333 


2.138 


366.502 


3.996 



Table 4. Same as Table |3] for the constant mass- to-light ratio models in lde Blok &: Bosmal 1 12002^ . 
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Figure 10. Continued. 
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responding best-fitting values of V200 and C200 (Table 4 in 
Ide Blok fc Bosmal 120021 ). Values of ps and t-b for the maxi- 
mum disk and constant mass-to-light ratio cases are given 
in Tables [3] and [4l respectively. 

Figures [9l and [TOl show the predicted rotation curves for 
70 = 200 (mdmC^ ~ 100 MeV) and 70 = 2000 (mdmC^ ~ 
1 GeV) with {crv)^± = 2.6 x 10"^^ cm^ s"^ and all other 
parameters set to our canonical values. Since these rotation 
curves have been computed (at least, in the innermost re- 
gions) from the Ha line, we also plot the results obtained 
for a gas density n = pg/nip = 0.01 cm~'^ and a;ion = 1, 
appropriate for the hot, diffuse compon ent responsible for 
the em ission line. The original model of Ide Blok fc Bosmal 
l|2002r ) without dark matter annihilation is shown for the 
sake of comparison, and the reduced x^ values associated to 
each model are listed in Tables [3] and [J] 

In general, the effect is not very significant for n = 
1 cm~^ (not even for rridmC^ ~ 100 MeV, except for a few 
exceptional cases, such as NGC 3274) . For the adopted value 
of the logarithmic slope of the dark matter density profile 
near the centre, a — 1, and the extre mely low values of 
the ch aracteristic density ps reported bv lde Blok fc Bosmal 
l|2002f) . the circular velocity at the innermost point becomes 
reduced by an amount that is typically much smaller than 
the observational error bars. A more noticeable effect would 
be obtained for steeper profiles (see Figure |8]), but also if 
one considers the typical density of the hot, ionized medium 
where the Ha line originates. Using n = 0.01 cm~^, the 
rotation curves of all galaxies would be dramatically affected 
on ~ kpc scales for a = 1 and Eo ^ GeV, both for the 
constant M/L and maximum disk models. 

These results represent a double-edged sword for dark 
matter annihilation models. On the one hand, it might be 
possible to find a particular dark matter candidate that is 
able to explain the rotation curve data with a cuspy den- 
sity profile. On the other hand, we also predict that, in that 
case, one should observe prominent differences in the kine- 
matics of the stellar, neutral, and ionized components due to 
their different densities. The observed rotation curves pro- 
vide thus an additional tool (complementary to radio and 
gamma-ray constraints) to rule out a broad class of models 
and hopefully help to identify the physical properties of dark 
matter particles. 



ter (intensity of the magnetic field, density and ionization 
fraction of the ISM gas, inner logarithmic slope of the dark 
matter density profile and the virial mass of the galaxy) to 
verify that our results hold in the general case. Our main 
conclusions can be summarized as follows: 

(i) For the canonical Milky Way model, the dark matter 
pressure gradient is able to offset the gravitational accelera- 
tion within the central ~ 10 — 400 pc as long as the injection 
energy is lower than 1 GeV. There would be an extremely 
weak signature if Eq ~ GeV, and the effect would be com- 
pletely negligible for larger values of -Eo- 

(ii) The ionization fraction of the ISM and the intensity 
of the magnetic field determine the energy losses and the 
shape of the electron-positron spectrum at low and high 
values of the Lorentz factor, respectively. Although these 
details may have a strong impact on other observables, such 
as the emission at different wavelengths, they do not affect 
the rotation curve significantly. The precise value of the gas 
density plays a more important role, and it changes the re- 
sults at the quantitative level. 

(iii) Steep logarithmic slopes of the dark matter density 
profile yield much higher pressures in the central regions. For 
Q J5 1.9, a clear signature of dark matter annihilation on the 
observed rotation curve is expected even for _Eo ~ 1 TeV. 

(iv) Comparison with publicly-available observational 
data shows that, while dark matter pressure may bring the 
predicted rotation curves into better agreement with obser- 
vations, it is arguably more likely that this effect is more 
useful as a constraint on the annihilation cross-section as a 
function of dark matter particle mass. 
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5 CONCLUSIONS 

In this paper, we have investigated the contribution to the 
total gas pressure arising from relativistic electrons and 
positrons produced in dark matter annihilations. The prop- 
agation of these particles through the ISM is determined 
by the diffusion-loss equation. We assume a uniform dif- 
fusion coefficient and consider inverse Compton scattering, 
synchrotron radiation. Coulomb collisions, bremsstrahlung 
and ionization of neutral hydrogen atoms as the main en- 
ergy loss mechanisms. All the electrons and positrons are in- 
jected with an initial energy Eo between 1 MeV and 1 TeV, 
and the injection rate is constrained by different Galactic 
observations. 

We have evaluated the effect of this "dark matter 
pressure" for astrophysical conditions representative of the 
Milky Way and varied the adopted values of each parame- 
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